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Abstract 

We compute the critical temperature T c of a weakly interacting uniform 
Bose gas in the canonical ensemble, extending the criterion of condensation 
provided by the counting statistics for the uniform ideal gas. Using ordinary 
perturbation theory, we find in first order (T c — T®)/T® = — 0.93a/} 1 / 3 , where 
is the transition temperature of the corresponding ideal Bose gas, a is the 
scattering length, and p is the particle number density. 

PACS numbers: 03.75.Fi 
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I. INTRODUCTION 



It may well look like a long-solved text book excercise, but the variation of the dilute 
Bose gas' critical temperature with the interaction strength has not yet found a conclusive 
answer. To date, all authors assume a continuous behavior in the limit of weak interaction, 
lim a ^ T c = T c °, where T c ° is the transition temperature of the non-interacting system, and 
a is the s-wave scattering length. However, the sign, the proportionality constant c, and the 
exponent 77 in the expression for the shift in the critical temperature at fixed density p 
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are still subject to considerable debate. Early calculations by Fetter and Walecka |l| and 
Toyoda 0] predict a decrease in temperature, T c < T c ° (Yet, it should be noted that the 
expression derived by Fetter and Walecka yields zero for a point potential.). However, more 
recent calculations indicate the opposite [3-10]. Concerning the exponent, one finds in the 
literature a set of predicted rational values which range from 77 = 1/6 ][| , ][J , [[U]] to r\ = 1/2 
1 12]] . The most recent analytical investigations converge towards the value t/ = 1/3 |], 0, 
||, 0, i. e. they predict a linear dependence of the temperature shift on the scattering 
length. This result is also backed by Monte Carlo simulations ||, ||, and by an ingenious 
extrapolation of experimental data for the strongly interacting condensed He4 in vycor glass 
|14| . Still, the result of Toyoda r] = 1/6 continues to find support ||. The proportionality 



constant, finally, has been predicted to assume a variety of values, which for rj = 1/3 range 
from c = 0.3 [§] to c = 5.1 ||14|| . The Paris group most recent numerical analysis, for 
example, points at c = 2.3 f9j, a value which is close to the theoretical prediction of Baym 
and collaborators [0, while the extrapolation of the experimental data on He4 in vycor glass 
favors c = 5.1, which is closer to an early prediction c = 4.66 of Stoof 

It is frequently maintained that ordinary perturbation theory can not be applied as it is 
plagued by seemingly unsurmountable infrared divergencies (See 0). We point out that this 
conclusion is based on the implicit assumption that the grand-canonical statistics, which is 
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governed by a chemical potential, is a sensible approximation to the real system, i. e. a system 
where - as a matter of principle - not the chemical potential, but rather the total number of 
particles is fixed, possibly at very large a value. While this assumption of thermodynamic 
equivalence does indeed hold in a system with sufficiently strong interactions, it must be 
rejected in the in the limit a — > 0. In this limit the grand-canonical statistics implies 
fluctuations of the ground state occupation, which for temperatures at and below T c turn 
out to be extravagantly large, An ~ 0(N) [15-17]. It is these unphysical fluctuations 
which doom to failure any attempt to reliably compute the shift of the Bose gas critical 
temperature, in the non-interacting limit a — > 0, when resorting to ordinary perturbation 
theory in the grand-canonical ensemble. 

As the ground state giant fluctuations are easily traced back to the fluctuations in the 
total number of particles, which in the grand-canonical statistics turn into an unacceptable 
AN ~ iV for T < T c °, a safe way out is to resort to statistical ensembles where the total 
number of particles is not allowed to fluctuate. In the canonical and microcanonical en- 
sembles, for example, the ground state fluctuations of the non-interacting system exhibit a 
scaling An ~ 0(N 2 ^) [15] which - although still anomalous - turns out to be sufficiently 
suppressed for ordinary perturbation theory to be applicable. 

Indeed, as we shall demonstrate in this letter, first order perturbation theory in the 
canonical ensemble yields the following shift in the critical temperature (where A is the De 
Broglie thermal wave length at T = T c °): 

8tt 
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which - contrary to some early expectations - is neither zero nor infinite. 
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II. THE HAMILTONIAN 



We consider a uniform system of N weakly interacting bosons in a volume V = L 3 , 
imposing periodic boundary conditions. The Hamiltonian reads 

H = H + H m t , (4) 

where H is the Bose gas kinetic energy, 

h 2 k 2 

Ho = V] e k h k , e k = — — , (5) 

2m 



k 

and H int describes the particle pair interaction, 



tt u ft ft i i Anh 2 aN 

pkq 

Here k = (2n/L)n is a wave vector, with n a vector of integers, 6k, are bosonic particle 
annihilation and creation operators, = is the associated number operator, m denotes 
the particle mass, and a denotes the s-wave scattering length. 



III. THE COUNTING STATISTICS 

We shall be working at fixed density p = N/V (equivalently: fixed specific volume v = 
p _1 ), but variable total number of particles iV (and, concomitantly, variable system volume 
V). The first issue to be faced is to provide a definition of the transition temperature, which 
- as we recall - only acquires the meaning of a critical temperature in the thermodynamic 
limit iV — > oo, V — > oo, p = N/V constant. 

We base our definition on the counting statistics of the zero-momentum state, 

Pn(P; N) = ^^yTr {6^^%^} , (7) 

which is the probability to find n particles (out of N total particles) in the zero-momentum 
state p = hk = 0. Here, N = J2k ^k is the operator for the total number of particles, 5 a ,b is 
the Kronecker delta, and Z(/3;N) is the canonical partition function, 
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Z(P;N) = Tr{e-P% iir } . (8) 
In the non-interacting limit, the counting statistics for high temperatures is a strictly 



decreasing function of n, i. e. P n > P n +i [18]. For sufficiently low temperatures it displays 



a single peak at n ~ (n) ~ O(N). Assuming that a system of weakly interacting bosons 
behaves correspondingly, we introduce the auxiliary function 

D(J3; N) = Tr {[<^, - S^i] ^h, N ) ■ (9) 

The cross-over from the high-temperature regime, where D > 0, to the low-temperature 
regime, where D < 0, is assumed to occure for a certain value /3 = /?*, which is defined by 
the relation 

D(J3.;N) = 0. (10) 

For fixed density p, and fixed scattering length a, the solution of this equation depends on 
the total number of particles N: /3* = fl*(N). We stipulate that, in the thermodynamic 
limit, the cross-over temperature T* = (^b/5*)" 1 coincides with the critical temperature of 
Bose-Einstein condensation: 

lim T*(N) =T C . (11) 

This identification, being non-trivial for an interacting system, will be verified below for the 
non-interacting case. 

IV. PERTURB ATIVE ANALYSIS OF T* 

We determine /3* using a series expansion in Hi n t- The Dyson series of D = D(/3; N) 
reads 

D = D + D 1 + D 2 + ... , (12) 
where _D n = D n (/3; N) is of n-th order in H m t. The first two terms are given by 



D (f3; N) = Tr {[5 hofi - S^e'^S^} 



(13) 



(14) 



To solve Eq. ( ]TD| ) we set 



& = + Aft 



(15) 



where f3^ denotes the cross-over inverse temperature of the non-interacting Bose gas, and 
A/3* is a correction which is assumed to be small. The defining equation for (3^ reads 



D O (0;N) = O 



(16) 



and the shift to leading order 



A/3* Di(P;N) 



PEa(P;N) 



P=P 



(0) 



(17) 



where 



EoiP; N) = ~D (P; N) . 



(18) 



A. Exact Relations 



Observing e — 0, which implies that Hq does not depend on no, we may recast Eq. ( |T3"D 
into the form 



D (P; N) = Tr ex { 



$N C ^N ^N ex ,N-l 



-PHo 



(19) 



where Tr ex denotes the trace over the occupation of excited states k ^ 0, and 



(20) 



denotes the operator of the number of particles in the excited states. Furthermore, using 

(21) 



N ! N Y 2 
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where R has no diagonal elements in the Fock basis, and observing that [Sn ,o — Sn ,i\n (n — 
1) = 0, one finds 

n k (n k - l) c _ H o 



uB 

A(AiV) = ^Tr cx 



(N-l)u(3D ((3;N). 

(22) 



Note that due to the definition of [3^ the second terms does not contribute to the shift A/?*. 
To proceed, we use the Laplace representation of the Kronecker delta 

1 r+iir 



S~ — — I dae( N ~ Ne * )a 
°n cx ,n ~ 2-jri 



(23) 



and perform the trace Tr ex - We then face 

_ 1 fin 

D ((3;N) = — da 



l-e- a 



,F(a) 



(24) 



DM N) = ^r^- jT da [l - e~ a ] I ]T n£(a)J — (N — l)u[3D ([3; N) , (25) 



AT 2ni j-t-K 
where n^a) = nk(«; /3, N), 

n k (a;[3,N) = -^— -, 

and F(a) = F{a;(3,N), 

F(a; [3, N) = Na + N^g 5/2 (a) . 
Here A = A(/3) is the thermal De Broglie wave length, 



(26) 



(27) 



X((3) = pnh 2 /(mk B T) , 

and (75/2(0:) = #5/2(0:; P, N) is a discrete predecessor of a Bose integral function, 

A 3 



(28) 



g 5/2 (a;(3,N) 



Nv 



(29) 



Upon identifying a = —flfJ,, where /x denotes the chemical potential in the grand-canonical 
ensemble, we note that for fixed a the corresponding value F(a) is nothing but the ideal 
Bose gas grand-canonical free energy, and n^{a) is the grand-canonical mean occupation. 



B. Continuum Approximation 



For sufficiently large N, and for the interesting range of thermal De Broglie wavelength 
such that A 3 p <C N, we may invoke the continuum approximation, and replace the discrete 



sum over momenta by an integral, J2k^o ~~ > J d 3 k.f\ We then obtain 



F(a; (3, N) — > F(a; (3, N) = Na + N—g 5/2 (a) 



]T K(a; (3, N) — > N— (g 1/2 (a) - 93/2(0)) 



where the g a (a) denote the Bose-Einstein integral functions 



g a (a) 



1 



dx- 



x 



a-l 



r(cr) Jo e x+a - 1 
Concomitant with the above replacements we also have: 

Do — D = [ £L \l - e -l e F ^ , 
PE -^f3E = N~h /2 (f3 ] N), 
Di — A = ^f3u [l 1/2 (/3; iV) - J 3/2 (/3; iV) 



(N-l)PuD . 



where 



W;N)= — 



da r 
2?ri 



fl'a(«)e 



F(a;/3,N) 



(30) 
(31) 



(32) 



(33) 
(34) 
(35) 



(36) 



C. Expansion in JV x / 3 



As A" is assumed to be large, N 3> 1, it is now tempting to evaluate Dq, Di in the saddle 
point approximation. Yet, due to the 0(N~ 2 ^ 3 ) proximity of the saddle point to the branch 



1 In an extended version of this paper we shall include a systematic study of the corrections to the 
continuum approximation. 
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point of F, this procedure is doomed to fail, and a completely different treatment must be 
developed.^ 

Since we expect the integrals to be dominated by the small values of a, we resort to the 



Robinson representation [19 



g a (a) = r(l - ^R- 1 + £ ^((a - n)a n . (37) 

n=0 n - 

Save for the branch cut of a a ~ l , which runs along the negative real axis, the Robinson 
expansion converges absolutely for \a\ < 2tt. Note that the radius of convergence covers the 
domain of integration in the above a-integrals. 

Exploiting the Robinson representation, the exponent reads 

F = In C - Ya + \xo? t2 + XO(a 2 ) , (38) 

where InC = [C(5/2)/(A 3 p)]iV is a constant, 0(a 2 ) denotes some analytic function, which 
may be extracted from Eq. (0), and 



V V? K N y 



X 3 

A 3 



N, (39) 



with Aq = C(3/2)/p, the thermal De Broglie wave length of the non-interacting gas evaluated 
at the transition temperature. 

We note that for given N, the solution of Eq. ([TDD implies a relation between X, which 
is O(N), and Y. As we expect A* — > Ao in the limit N — > oo, the scaling of Y* with N is not 
obvious, yet 

e=J (40) 
will certainly be small. Introducing the transformation of the integration variable, 

a — > r = e _2 o; (41) 



2 In field theory, the proximity turns into a confluence which renders meaningless the non- 
interacting limit of the theory at T = Tq. 

9 



the free energy reads 



F = lnC + A 



where 



" T+ 3 r3/2 ) + Aer 5/2^; e ) ; 



A Y " 



(42) 



(43) 



and r 5 / 2 is a regular function, 

^5/2 (^;e 



r 2 ^ (-e 2 Y 
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C(i/2-^K 



(44) 



Since we shall find that A^ ~ 0(1) at the cross-over temperature T = T*°\ and 
concomitantly el°^ ~ 0(N^ 1 ^ 3 ) for large AT, we may invoke a formal expansion in e, 



L) " e 4 ^(A) + e 5 ^^AK 3 (A) + 0(e 6 ) , 



C 



40F 



(45) 



71/2 " 13/2 " C(V2) " C(3/ -^^o + e^K 1/2 (A) + e^A^ /2 (A) + 0(e 5 ) , (46) 



C 



C 



^ = ^A, - e\(3/2)K 2 (A) - J 



C(3/2)C(l/2) . v 4v^F 



40F 



AAT 4 - 



5/2 



+ 0(e 8 ) , (47) 



where we have introduced the family of functions 



(48) 



The functions K v obey the recurrence relation 



Kv+3/2 = K v 



+i 



A 



K v 



(49) 



which is easily proven by expressing K in terms of X and Y, using the inverse of the 
transformation fl4lT). 
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D. Results 



Upon inserting Eq. into Eq. flTfll) , the condition which fixes the cross-over temper- 
ature of the non-interacting gas reads 

^i(Al 0) )=0. (50) 

up to corrections 0(N^ 1 ^ 3 ). This equation is easily solved numerically, yielding the result 
a!°^ = 0.334. Expressed in terms of temperature we have 



T (o) _ T o 



(o) \ V3 



( 32ttA1 
+ 1 27C(3/2) 2 I A" 1 :! 



C(iV- 2 / 3 ) 



(51) 



where T c ° is the critical temperature of the non-interacting Bose gas. Note that for N finite, 
the cross-over temperature is slightly higher than the transition temperature of the ideal 
Bose gas, but in the thermodynamic limit N —>■ oo they coincide. 

Collecting terms and observing that D (/3i°^; N) = 0, the interaction induced shift reads 



Ap, 8tt 



^ 1/2 (A) 



ja=a!°) 



(52) 



0® 3C(3/2) A (0) [AK 2 (A) 

up to corrections of order 0(N~ 1 ^ 3 ). The shift involves the ratio / = Ki/ 2 / (AK 2 )\a=a* 
Exploiting the recurrence relation ([S), and observing that Ki(A*) = 0, we find / = —2/5 
We note that this value is exact as it does not depend on the numerical value of A* . 



We are now in the position to consider the thermodynamic limit of Eq. (|52j) . Identifying 
limjv^oo —T c , the result reads 

AT C _ T C - T e ° 2 8tt a 

T° _ TO 5 3C(3/2)A 1 ' 

-lmkr^ m -°^ 1/3 ■ < 54 » 

We thus find a negative shift in the critical temperature, growing linearly with the 
scattering length. This result can be compared with the other fully analytical prediction 
existing in the literature, derived by Baym, Blaizot and Zinn- Justin [0: 

AT C 8tt 



T c ° 3C(3/2) 4 /3 
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ap 1/3 = 2.33ap 1/3 . (55) 



The two results exhibit the same scaling, but differ for the sign and for a factor 2/5 in 
the proportionality constant. 

The prediction of Baym et al. has been obtained by evaluating the leading order in 
the 1/N expansion for a O(N) field theory model that coincides with the original Bose 
Hamiltonian for N = 2, and by observing that the final result does not explicitely depend 
on N. However, the result is strictly proven only for large N, and whether it is reliable also 
for N = 2 is still an open problem. 

On the other hand, our approach based on the counting statistics and on ordinary per- 
turbation theory in the canonical ensemble indicates that in the limit of ultraweak inter- 
action there are contributions, otherwise possibly neglected in other approaches, that tend 
to suppress quantum effects. Clarification of this issue may be expected from higher order 
perturbation theory. 
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